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ABSTRACT 

We have modeled the time- variable profiles of the Ha emission line from the non-axisymmctric disk and 
debris tail created in the tidal disruption of a solar-type star by a 10^ black hole. Two tidal disruption 
event simulations were carried out using a three dimensional relativistic smooth-particle hydrodynamic 
code, to describe the early evolution of the debris during the first fifty to ninety days. We have calculated 
the physical conditions and radiative processes in the debris using the photoionization code CLOUDY. 
We model the emission line profiles in the period immediately after the accretion rate onto the black 
hole became significant. We find that the line profiles at these very early stages of the evolution of the 
post-disruption debris do not resemble the double peaked profiles expected from a rotating disk since 
the debris has not yet settled into such a stable structure. As a result of the uneven distribution of 
the debris and the existence of a "tidal tail" (the stream of returning debris), the line profiles depend 
sensitively on the orientation of the tail relative to the line of sight. Moreover, the predicted line profiles 
vary on fairly short time scales (of order hours to days). Given the accretion rate onto the black hole we 
also model the Ha light curve from the debris and the evolution of the Ha line profiles in time. 

Subject headings: black hole physics — galaxies: nuclei — line: profiles 



1. INTRODUCTION 

1.1. Tidal Disruption of a Star by a Black Hole and 
Related Issues 

A star in an orbit around a massive black hole can get 
tidally disrupted during its close passage by the black hole. 
After several orbital periods the debris from the disrupted 
star settles into an accretion disk and gradually falls into 
the black hole (Rees 1988; Cannizzo, Lee & Goodman 
1990; Syer & Clarke 1992; Loeb & Ulmer 1997). As mate- 
rial gets swallowed by the black hole intense UV or soft-X 
ray radiation is expected to emerge from the innermost 
rings of the accretion disk (Frank & Rees 1976; Lightman 
& Shapiro 1977; Frank 1978; Phinney 1989; Sembay & 
West 1993; Magorrian & Tremaine 1999; Syer & Ulmer 
1999). For black hole masses Mbh < 1O''M0, tidal disrup- 
tion theory predicts flares with luminosities of the order 
of the Eddington luminosity with durations of the order 
of months, and with spectra that peak in the UV/X-ray 
domain band (Rees 1988; Evans & Kochanek 1989; Ulmer 
1999; Kim et al. 1999; Gezari et al. 2002). High-energy 
flares from the central source illuminate the debris, the 
photons get absorbed, and some are re-emitted in the op- 
tical part of the spectrum (i.e. the light is "reprocessed"). 
One of the spectral lines in which this phenomenon can be 
observed is the Balmer series Ha line {Xrest = 6563 A). 

The disruption of a star begins when the star approaches 
the tidal radius, ~ r^,(M;,;i/Mj,)^/'^, the point where the 
surface gravity of the star equals the tidal acceleration 
from the black hole across the diameter of the star (r* 
and M^, are the radius and mass of the star and Mth is 
the mass of the black hole). A 10^ Mq black hole is often 
used as a prototypical example in tidal disruption calcu- 
lations. This choice is motivated by the criterion for a 



solar-type star to be disrupted before it crosses the black 
hole horizon (i.e. the Schwarzchild radius, r^) in order for 
emission to be observable. For supermassive black holes 
with Mbh > 10® Mq, rg > rt and the star falls into the 
black hole before it gets disrupted. 

The tidal disruption process has been the subject of 
many simulations (Carter & Luminet 1982, 1983; Bick- 
nell & Gingold 1983; Evans & Kochanek 1989; Khokhlov, 
Novikov, & Pethick 1993; Laguna et al. 1993b; Frolov et 
al. 1994; Marck, Lioure, & Bonazzola 1996; Deiner et al. 
1997; Ayal, Livio, & Piran 2000; Ivanov & Novikov 2001; 
Ivanov, Chernyakova, & Novikov 2003). It has been shown 
that tidal processes in the vicinity of a massive black hole 
could lead to tidal capture, tidal heating and tidal spin-up 
of a star (Novikov, Pethick, & Polnarev 1992; Alexander 
& Kumar 2001; Alexander & Hopman 2003; Alexander & 
Morris 2003) , and in some cases ultimately to the explosion 
of the star. Such explosions, as well as accretion of post- 
disruption debris, should manifest themselves as luminous 
flares from the centers of galaxies (Carter & Luminet 1982; 
Rees 1988). Stars close to a black hole may experience 
mixing or may eject some of their mass (Alexander 2001). 
As a consequence, stellar populations in nuclear clusters 
are expected to be somewhat unusual in comparison with 
populations whose evolution was not affected by a mas- 
sive black hole (Alexander & Livio 2001; Di Stefano et 
al. 2001). This has important implications for observa- 
tions of stellar cluster in the center of our Galaxy (Chez 
et al. 2000; Schoedl et al. 2002; Eckart, Ott, & Genzel 
1999; Figer et al. 2000; Gezari et al. 2003) where high con- 
centration of otherwise rare blue He supergiants has been 
observed (Krabbe et al. 1991; Najarro et al. 1994). 

The tidal disruption and accretion of stars can fuel 
black holes in the centers of galaxies (Hills 1975; Dun- 
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can & Shapiro 1982; David, Durisen, & Cohn 1987a,b; 
Murphy, Cohn, & Durisen 1991; Freitag & Benz 2002; Yu 
2003) and its contribution to nuclear activity in galax- 
ies and the growth of the black hole mass depends on 
the rate of disruption events in a galaxy. The predicted 
tidal disruption rate in a typical inactive galaxy is 10~^- 
IQ-^ yr"^ (Magorrian & Tremaine 1999; Alexander 2001). 
This value is consistent with the rate of UV/X-ray out- 
bursts observed with RO SAT from inactive nuclei selected 
as tidal disruption candidates (Donley ct al. 2002). The 
rate of tidal disruptions in active and more luminous nu- 
clei is estimated to be lower, with the lowest value of 
10"''' yr^^, for galaxies with the most massive black holes. 
This may occur partly because massive central black holes 
{Mbh > 10^ Mq) swallow stars promptly, without disrup- 
tion, and partly because stars are less centrally concen- 
trated in these galaxies (Magorrian fc Tremaine 1999). 
The observed UV/X-ray flaring rate in these galaxies is 
about 9 X 10~^ galaxy"^ yr~^ (Donley et al. 2002) suggest- 
ing that in such nuclei, outbursts may be due to another 
mechanism, such as accretion-disk instabilities (Sicmigi- 
nowska, Czerny, & Kostyunin 1996; Burderi, King, & 
Szuszkiewicz 1998). 

The tidal encounter of a compact star with a black hole 
can also result in emission of gravitational waves, which 
may be observable with upcoming instruments. More 
specifically, compact stars (helium stars, white dwarfs, 
neutron stars, and stellar-mass black holes) which can 
withstand large tidal forces without being disrupted, may 
get captured in relativistic orbits around a supermassive 
black hole. Due to the in-spiral and decay of the orbit 
those objects arc expected to emit the peak of their grav- 
itational wave power in the LISA frequency band (Hils & 
Bender 1995; Sigurdsson & Rees 1997; Freitag 2001). It 
has been recently suggested by Freitag (2003) that very- 
low mass main sequence stars (MSSs; Af <C 1 Mq) may 
contribute to events detected by LISA, which was not pre- 
viously expected for capture of these objects by a super- 
massive black hole. Although MSSs produce a relatively 
weak gravitational signal during the in-spiral, compared 
to compact objects, their detection in the Galactic center 
is more likely because such stars have a predicted close- 
encounter rate that is an order of magnitude higher than 
that of white dwarfs (WDs) , and about two or more orders 
of magnitude higher than that of neutron stars (NSs) and 
stellar- mass black holes (BHs). These compact MSSs are 
expected to produce a strong enough signal to allow for 
0.5-2 detections from our Galactic center, with a signal- 
to-noisc ratio of 10 or higher, for a LISA mission duration 
of one year (Freitag 2003). Moreover, MSSs are expected 
to be disrupted relatively early during the in-spiral, giv- 
ing rise to possibly detectable electromagnetic flares. The 
sudden appearance of an electromagnetic counterpart to a 
transient gravitational wave source, expected in the case 
of MSSs and helium stars, could allow identification of the 
progenitor. In the case of a tidal disruption event in an- 
other galaxy, the coincidence of an electromagnetic flare 
and a gravitational wave signal would provide an indica- 
tion that the event occurred at the very nucleus of the 
galaxy, and possibly allow the measurement of the red- 
shift. More compact objects which can spiral in without 
being disrupted (such as WDs, NSs, and BHs) are ex- 



pected to produce stronger gravitational wave signatures 
and weaker or no electromagnetic flares, with the excep- 
tion of white dwarfs in which tidal interaction may trigger 
thermonuclear explosion (Garcia-Senz, Bravo & Woosley 
1999). 

1.2. Observational Motivation: Transient Emission 
Lines in Inactive Galaxies and LINERs 

In view of the above theoretical considerations, it is nec- 
essary to make predictions of the likely observational sig- 
natures of a tidal disruption event. The prompt UV/soft- 
X-ray flash that is expected to accompany the disruption 
does provide strong evidence for such an event and has in 
fact been detected in a number of galaxies with ROSAT 
and HST (Brandt, Pounds & Fink 1995; Grupe et al. 
1995a, b; Donley ct al. 2002; Bade, Komossa. & Dahlcm 
1996; Komossa & Greiner 1999; Greiner et al. 2000; Grupe, 
Thomas, & Leighly 1999; Renzini et al. 1995; Li, Narayan, 
& Menou 2002). However, the duration of this flash is 
short enough that it can easily be missed. Aftereffects with 
a longer duration, such as line emission from the debris, 
have a better chance of being detected. If the appear- 
ance of emission lines just after an X-ray flare could be 
detected from the same object it could be used to identify 
the tidal disruption in the early phase and would provide 
strong support for the overall picture, but such cases are 
rare (see Cappellari et al. 1999; Gezari et al. 2003). 

A set of tantalizing observations in the past decade 
show that several LINERs (low-ionization nuclear emis- 
sion regions; Heckman 1980) have transient Balmer emis- 
sion lines, which are often double-peaked; examples in- 
clude NGC 1097 (Storchi-Bergmann, Baldwin & Wilson 
1993), M81 (Bower 2000), NGC 4450 (Ho et al. 2000), 
NGC 4203 (Shields et al. 2000) and NGC 3065 (Eracleous 
& Halpern 2001). Such line profiles are characteristic of 
rotating disks and resemble the persistent double-peaked 
Balmer lines found in about 10-20% of broad-line radio 
galaxies and about 3% of all active galaxies (e.g., Era- 
cleous & Halpern 1994, 2003; Stratcva et al. 2003). Their 
abrupt appearance in LINERs led to suggestions that this 
transient event was related to the tidal disruption of a 
star by a supermassive, nuclear black hole (Eracleous et 
al. 1995; Storchi-Bergmann et al. 1995) or a change in the 
structure of the accretion disk associated with a change in 
accretion rate (Storchi-Bergmann et al. 1997). 

To investigate the possibility of line emission from the 
post-disruption debris and to evaluate the suggestion that 
the transient double-peaked lines of LINERs are related to 
tidal disruption events, we have undertaken a calculation 
of the strength and profile of the Hq line emitted from the 
debris. In §2, We describe two SPH (Smoothed Particle 
Hydrodynamics) simulations of tidal disruption on which 
we base our further calculation of the line properties. Our 
line profile calculation follows the method used for line 
profiles emitted by relativistic, Keplerian disks and is de- 
scribed in §3. In §4 we present the resulting line profiles 
and in §5 we discuss the physical conditions in the debris 
and the approximations used. In §6 we summarize our 
conclusions and consider future prospects. 

2. QUALITATIVE DESCRIPTION OF TWO SPH 
SIMULATIONS 
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Tidal disruption simulations were carried out with a 
three dimensional relativistic SPH code in order to study 
the dynamical evolution of the post-disruption debris. The 
SPH code used provides a description of relativistic fluid 
flows in a static curved spacctimc geometry (Laguna ct al. 
1993a, b). We use it to simulate the tidal disruption of a 
star in the potential of a Schwarzchild black hole. 

The main sequence star is modeled as a polytrope with 
index F = 5/3. The density profile of the pre-disruption 
star is determined by the Lane-Emdcn equations. The 
star is initially placed on a parabolic orbit at a distance of 
700 Tg from the black hole where =rs/2 = GM/c^ = M 
is the gravitational radius and M is the mass of the black 
hole. Hereafter we use units in which G = c = 1, where 
G is the gravitational constant and c is the speed of the 
light and we adopt = M as a natural unit of length. 
The gravitational radius of a 10^ black hole is Vg = 
1.48 X 10^-'^ cm = 4.92 light seconds and the dynamical time 
at a given radius is Tdyn ^ {r^/GMy/^ = 4.92 (r/rg)^/^ s. 

The strength of the tidal encounter, is given by the ratio 
of tidal radius to pericentric distance, r]=rt/rp . For the 
case of 1 M0 star and a 10^ M0 black hole n ~ 47 r^. 
The two SPH simulations describe the tidal disruption for 
the case of mildly relativistic encounter, ?7=1.2. This value 
has been selected as a likely scenario for tidal disruption of 
a main sequence star. We do not investigate the strongly 
relativistic cases where tidal compression could lead to the 
explosion of a star, because such an explosion could intro- 
duce an uncertainty in the distribution of the debris mass 
over binding energy and consequently in the spatial distri- 
bution and kinematics of the debris. 

The self gravity of the star is accounted for initially. 
Once the star gets disrupted the self-gravity becomes 
unimportant and the debris particles follow nearly Kep- 
lerian orbits. It can be shown for the case of the Keple- 
rian potential that the rate of return of bound debris to 
the pericenter follows dM/dt oc t~^^^ (Rees 1988; Phinney 
1989). This behavior of the debris return rate has been ob- 
served in our SPH simulations. Once bound debris starts 
to rain down on the black hole it is expected to cause the 
initial rapid rise in the emitted UV/X-ray light curve and 
steady decay with the power law index of — % later on. 

The two different simulations have 5,000 and 20,000 
particles (hereafter 5k and 20k respectively) contributing 
equally to the mass of a 1 Mq star. The 5k SPH calcula- 
tion follows the debris for 94 days in total. After 34 days 
significant accretion onto the black hole begins. Our in- 
vestigation follows the evolution of the line profiles in the 
last 60 days. The 20k SPH simulation spans 53 days, dur- 
ing which the evolution of the line profiles was followed for 
the last 6 days (Table 1). Using both the 5k and 20k SPH 
simulations in the line profile modeling we take advantage 
of the longer time span in the former and better resolution 
achieved with the larger number of particles in the latter. 

Figure 1 shows particle distribution maps after the sec- 
ond pericentric passage, at the beginning of the accretion 
phase and at the end of the 5k simulation. At the early 
stages of the tidal event most of the particles were located 
in the pronounced tidal tail. Sixty days later, about 20% 
of the particles are scattered from the tidal tail and form a 
quasi- spherical distribution, with most of its mass concen- 
trated in the equatorial plane (Cannizzo et al. 1990; Loeb 



& Ulmer 1997; Ulmer, Paczynski & Goodman 1998; Ulmer 
1999; Menou & Quataert 2001). This is a consequence of 
the intersection of the leading part of the tidal stream with 
itself (Kochanek 1994; Lee, Kang, & Ryu 1996; Kim et al. 
1999; Ayal ct al. 2000). We refer to the spheroidal part of 
the debris as the halo and to its planar component as the 
disk. The remaining 79% of particles are still confined to 
the tail and 1% are accreted onto the black hole. There 
is a concern that some fraction of particles (^ ^^VlV, 

where and N are the total duration and total number 
of particles in the run) of the halo are an artifact of the 
SPH simulation. This can be due to the tendency of the 
SPH numerical method to preserve the constant number 
of neighbor particles for each particle during the calcu- 
lation. In the regions with a small density of particles 
this leads to a "smoothing" over a large spatial range and 
it may introduce the scatter of particles from the debris 
plane to the halo. These particles cannot be distinguished 
from the population of particles scattered out of the debris 
plane by the intersection of the tidal tail with itself. We 
further discuss the implications of the spheroidal halo for 
the emission line profiles and total Ho; huninosity in §5.2. 

The velocity distribution in the tail is "bimodal" where 
the central part of the tail exhibits very low radial veloc- 
ities: particles on the near (right-hand) side of the tail 
"flow" towards the black hole, while particles on far (left- 
hand) side are moving in the opposite direction. This is a 
consequence of the energy distribution throughout the de- 
bris in the disruption process: after the initial disruption 
event 50% of the debris stays bound to the system and 
the other 50% is unbound. This effect has been predicted 
by theory in the case of stellar disruption after a single 
fly-by of the star (Rees 1988) and has been observed in 
tidal disruption simulations. The symmetry in the distri- 
bution of the debris over binding energy is a consequence 
of the spin-up of the star at the expense of orbital kinetic 
energy. The spin-up initially causes the development of 
a quadrupolar deformation. As the tidal interaction gets 
stronger the star starts to shed its mass since the material 
in the stellar bulge has reached the escape velocity at the 
star's surface. One portion of the stellar debris ends up 
deeper in the potential well of the black hole, which causes 
further spread in binding energy of the debris. This effect 
determines which portion of the debris stays bound to the 
black hole Rees (1988). Following the second passage of 
the debris through pericenter approximately 66% of the 
mass is unbound, 33% remains bound and only about 1% 
is accreted by the black hole. The maximal approach- 
ing and receding velocities in the debris, with respect to 
the stationary observer positioned at infinity, are of order 
10~^c. The dynamical evolution in the 20k run is the same, 
it is just followed over a shorter evolutionary time scale. 
Since the tidal tail includes a large fraction of boimd and 
unbound particles in both simulations, its morphology and 
velocity field greatly influence the observed line profiles. 

3. LINE PROFILE CALCULATION AND TIME DELAY OF 
REPROCESSED LIGHT 

3.1. Calculation of Line Profiles 

We follow the line proflle calculations carried out by 
Chen & Halpern (1989) and Eracleous et al. (1995) to ob- 
tain the observed profile from a Keplerian, relativistic, thin 
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disk in the weak field approximation. The description of 
the debris as a flat, thin structure is justified by the fact 
that the height of the debris is three orders of magnitude 
less then its dimensions in the orbital plane. The main ob- 
jective of the calculation is to obtain the final expression 
for the flux density in the observer's frame as a function 
of parameters defined in the reference frame of the debris. 
Figure 2 shows the coordinate system and the geometry 
of the debris. The observer is located on the positive z- 
axis, at a distance d +00 and above the orbital plane at 
i = 30° to the z'-axis. Since the calculation is presented in 
cited papers we just introduce its main steps and comment 
on its application to the case of tidal disruption debris. 

The total emission-line flux received from the debris by 
an observer at infinity is given by an integral over the plane 
of the image produced at infinity, namely 



d^h 



(1) 



where v, Iv, and dO. are the frequency, specific intensity 
and solid angle element measured in the frame of the image 
(i.e., of the observer). Using the impact parameter of rays 
at infinity, 6, as a coordinate in the image plane and ex- 
ploiting the Lorentz invariance of the quantity Iv/v^ and 
the fact that the debris is confined to a plane, equation (1) 
can be transformed into an expression for the flux density 
(i.e., the line profile) in terms of coordinates and physical 
quantities in the source frame (see detailed derivation in 
Chen, Halpern, & Fihppenko 1989, Chen & Halpern 1989, 
and Eracleous et al. 1995) 



dfll^ 



M'^vq cosz 



t>' h^D\i,4>')i,{^,4>') (2) 



where the new polar coordinates in the debris plane are the 
dimensionless radius ^ = r /rg and the azimuthal angle <j)' . 
In practice, the integration is performed by summing over 
particles, assigning to each particle an emissivity according 
to its position, as derived from calculation with the code 
CLOUDY (Ferland 1996, see §4.1). The limits of integra- 
tion describe the portion of the debris that emits the Ha 
line, between radii and ^out- The function D[^,(f>'), 
the "Doppler factor," is determined by the phase space 
distribution of the emitting particles and the metric, and 
describes the effects of gravity and the motion of the emit- 
ting particles on the energies of the emitted photons. The 
function 'ip{^,(j)') is determined by the geometrical distri- 
bution of the emitting particles and the metric, and de- 
scribes the effects of curved trajectories of light rays. In 
the special case where the debris is confined to a plane, 
these functions are given by 

^ ^ (1-2/OV^ L_f3r'[l- {h/rf{l - 2/0 f 
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with 7 the Lorentz factor, /?,■' and f3^' the radial and az- 
imuthal velocities of debris particles in the source frame, 
and 



, ,^ ,,, , 1 /I — sin 2 cos( 

^ ^,<^' =1 + 7 L ^ . . — 
4 \ 1 -|- sm z cos ( 



(4) 



The above analytic expression for ?/'(^,(/)') has been de- 
rived in the weak field approximation and is accurate to 
order This approximation is appropriate in our case 
because the portions of the debris that experience a strong 
gravitational field are also highly ionized and make a neg- 
ligible contribution to the Ha fiux (see discussion in §5.1). 
The ratio b/r describes how rays emitted from the de- 
bris are mapped to points in the image at infinity and 
is given by b/r = (1 — sin^ icos^ Finally, 
the Lorentz factor is given by 7 = [1 — /3^/(l — 2/^)~^ — 

Pi,ii-2/o-'r'/'. 

The emission properties of the debris are described by 
the local specific intensity, We take the local line 

profile to be a Gaussian corresponding to a velocity dis- 
persion a (in units of the speed of light). The width of 
the local line profile represents not only internal motions 
of the line-emitting gas (not captured by the SPH simu- 
lation) but also the velocity range between the discretized 
debris points we use in our numerical integration. We fur- 
ther assume that the emissivity of the line is a power-law 
with radius (see discussion in §4.1 and §5.1). Therefore, 
we write the specific intensity as 



2 (27r)3/2 a 



exp 



{1 + X-Df 
2£>V2 



(5) 



where X is defined hy 1 + X = u/uq, where u and vq are 
observed and rest frame frequency, and eo is a constant. 

In the final model, the line profiles are described by the 
following parameters: the inner and outer radius of the 
line-emitting portion of the debris, and ^out, the parti- 
cle emissivity power-law index g*"*' and (jr***' , the inclina- 
tion of the debris plane z, the local velocity dispersion cr in 
units of c, and the radius of the central continuum source 
^0- The last parameter is relevant to the calculation of the 
light travel time across the debris, which we describe in 
§3.2, below. 



3.2. Time Delay of the Reprocessed Light Emitted by the 

Debris 

In the emission model we adopt, a central source of finite 
dimensions illuminates the debris. The luminosity of this 
source is proportional to the accretion rate onto the black 
hole. Due to the finite velocity of light, at any given time 
(in the observer's frame), difi^erent portions of the debris 
are seen to respond to a different level of illumination. It 
is also noteworthy that the length scales and time scales 
in this problem span a very large dynamic range. As a 
consequence, the light-crossing time of the outer portions 
of the debris is longer than the dynamical time of the inner 
portions, which makes it necessary for us to follow the re- 
distribution of the debris in phase space and the variations 
of the X-ray source carefully. 

It is simple to show that travel time delay for light rays 
can be written as 



^ttravel = (^ " ^o)(l - sin i COS^') 



(6) 



where ^ and (f>' are the coordinates of a particle in the 
orbital plane, ^0 is the radius of the central source and i 
is the inclination. Since light rays travel in the gravita- 
tional potential of the black hole they suffer an additional, 
relativistic time delay, which can be calculated from the 
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equation of geodesies for photons (Weinberg 1972, p. 202). 
The assumptions are that photons travel in an isotropic 
gravitational field and that their trajectories can be con- 
sidered coplanar with the observer and the black hole. We 
consider only the gravitational delay caused by the black 
hole's potential and assume that the debris does not have 
any significant gravitational influence on a light ray. This 
is a reasonable assumption since pdebris/Pbh ~ 10~^^. In 
our notation, the gravitational time delay for non-rotating 
black hole can be expressed as 



At 



gr 



1/2 



2 In 



1/2' 



^cosi 



- COS( 



1 — coszV^^ /l-|-sini 

— : + 2 In 

1 + COS t J \ COS I 



(7) 



The overall delay for a particle caused by light travel and 
general relativistic effects is then 



At = At 



travel 



Ato 



Ogr- (8) 

We only need to account for relative time delays, which we 

calculate relative to a ray coming from the origin of the 
coordinate system of the debris. We find that the general 
relativistic time delay in our calculations never exceeds 
10% of the travel time delay, and it is typically of order a 
few percent. 

4. RESULTS AND IMPLICATIONS 

4.1. Light Curves and Emissivity of the Debris 

The output of the 5k run comprises 351 frames describ- 
ing the evolution of the debris morphology over 60 days of 
accretion, with a time step of approximately 4 hours (see 
Table 1) . This is a fine enough temporal resolution to trace 
the redistribution of particles in the debris. For compari- 
son, a particle at r = 200 rg orbits the black hole with a 
period of approximately 24.5 hours. During the accretion 
phase of the simulation the number of particles decreases 
due to infall in the black hole. Tracking this number al- 
lows us to follow the accretion rate and construct the X-ray 
light curve. The total amount of mass accreted in this run 
is about 10"^ Mq. The X-rays resulting from accretion 
illuminate the debris out to large distances from the black 
hole and power the emission of Ha photons. 

The output of the 20k run consists of 204 temporal 
frames spanning 6 days of accretion with a time step of 
about 45 minutes (the coverage is not even; it includes 
gaps since the behavior of the debris can be captured even 
with sparse sampling). Because of the short time span of 
accretion, the number of accreted particles is small (less 
than 0.1% of the total mass) and consequently, the illu- 
minating, X-ray light curve is not smooth. To address 
this issue we compute the Ha light curve for two more, 
fiducial, illumination light curves: one consistent with the 
debris return rate predicted by theory (cx t"'''/''; Rccs 1988: 
Phinney 1989) ^, and the other, constant in time. This al- 
lows us to also examine the effect of different illumination 
patterns on the line profiles, which we describe in §4.2. 

In order to calculate the resulting Ha luminosity of the 
debris it is necessary to determine the efficiency with which 

^ We have verified that the energy distribution of the debris particles 
pericenter, based on the Keplerian orbits 



the debris reprocesses the incident radiation. This effi- 
ciency can be characterized by the surface emissivity of 
the debris as a function of radius, e = eo as shown in 
equation (5) . We used the photoionization code CLOUDY 
(Ferland 1996) for numerical calculations of the response 
of the debris to illumination. From the calculated sur- 
face emissivity as a function of radius (ergs~^cm~^) we 
were able to determine the particle emissivity as a function 
of radius Q = e/E = Qo C^'' (erg s~-^ particle"-^), where 
P = q — p is the particle emissivity power-law index, given 
the surface density distribution of the debris, S oc 

The response generally depends on a spectral energy dis- 
tribution (SED) of the incident radiation, matter density 
distribution, and the system geometry. We have described 
the debris system in terms of these physical parameters 
and performed CLOUDY calculations. 

We have assumed a SED for the illuminating radiation 
of the form oc i^^^ , extending to 100 keV. The nor- 
malization for the SED is set by the accretion luminos- 
ity of the debris. The mass accreted during the accretion 
phase of the 5k simulation gives rise to a time-average 
accretion luminosity of Lace ~ 1-5 x 10^^ Mq erg s~^, 
while the maximum luminosity achieved at the begin- 
ning of the accretion is L^^c^ « 8.0 x 10'^'* Mq erg s~^ 
(where Mq = Mfe/j/lO^M©). The luminosity during flares 
therefore may exceed the Eddington luminosity, LEdd = 
1.51 X 10^^ Me erg s'^. For the purpose of CLOUDY 
calculations, in order to charac;terize the emissivity of the 
debris, we have used the time-average value for the ac- 
cretion luminosity. Our estimated accretion luminosity is 
comparable to what was observed by Gezari et al. (2003) 
for three tidal disruption candidates, namely Lx = 10**^- 
10'"' erg s-i. 

Based on the particle distribution in the simulation we 
find that the density in the tidal tail decreases with time 
because the tail gets stretched as the debris evolves (cf.. 
Figure 1). We find that for the typical SPH time frame 
the density distribution in the tail can be approximated 
as n^*' oc ^ and is in the range n^*' = 10^^-10^^ cm~^. 
With path lengths of 9 x 10^^ cm and ~ 5 x 10^^ cm, 
the corresponding tail column densities are in the range 
j^taii ^ io27_io3i cin-2. The low and high values of the 
column density correspond to lines of sight along the short 
and long axes of the tail, respectively. (See Table 2 for 
values of physical parameters in the debris.) The density 
of particles scattered from the tail in the spheroidal halo 
around the black hole is fairly low in comparison. The 
halo density reaches a maximum in the plane of the de- 
bris (> 3(7 over-density), where particles orbiting around 
the black hole form a disk of radius ~ 2500 Vg. The num- 
ber density in the disk decreases with radius as nj}"'' = 
1 X 10^^ cm~^(^/500)~^-^. The corresponding column den- 
sities in the directions orthogonally and radially through 



the disk are Ar*^'= s:^ 7 x 10^° - 3 x 10^'"' cm^^, respectively. 
The spheroidal part of the halo, formed from particles scat- 
tered out from the plane of the debris has the estimated 



number density of n]f 



halo 



3 X 10 
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ne/500) 



-1.4 



and 



1 X 10^5 cm-2. 



corresponding column density of iV^"'° 

The CLOUDY calculations show that the Ha power 
emitted by the mostly neutral tidal tail in response to X- 

in our simulation would lead to a oc return rate of particles to the 
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ray illumination decays with distance from the center as 
gtaij ^ 7.-1-6^ ag a, consequence of the density distribution 
in the tail and the geometry of the debris (i.e. illumina- 
tion incidence). The corresponding particle emissivity for 
the tail derived from the surface emissivity decays with 
the distance as Q*"*' oc r~^-^ (Fig. 3). The disk and the 
halo components are almost completely ionized and only 
1 in 10^ hydrogen atoms is neutral, on average. The emis- 
sivity of the disk decreases with distance from the black 
hole as g''**'^ cx r^^'^, approximately. The equivalent par- 
ticle emissivity distribution for the disk is Q''**'^ oc r~'' °^. 
Therefore, in our description of the surface emissivity 
of the debris (equation 5) we use power-law indices of 
qtail _ 2 5 for the tail and g'^*"'^ = 2.2 for the disk. In our 
calculation for line profiles and light curves, due to the dis- 
crete nature of the SPH simulation, we use the prescription 
for the emissivity per particle, with power-law indices of 
ptaii ^ 2.4 for the tail and = 0.06 for the disk. The 

results of the CLOUDY calculations and the power-law 
prescriptions derived from them are summarized in Fig- 
ure 3. The total Ha luminosities contributed by the tail 
and the halo for the time-average value of the illumination 
Lace = 1.5 X 10'*'' erg s^-*, as calculated by CLOUDY, are 
^tc»( ]^o36 erg s^i and L*^'' ~ 1 x lO^"^ erg s'K The 

calculated value of the luminosity for the spheroidal com- 
ponent is ~ 6 X 10^* erg s^* (also hsted in Table 2). 

We discuss the emission properties of the debris in §5.1 
and their implications for the observability in §5.3. Fi- 
nally, with the above emissivity prescriptions we calculate 
the observed Ha luminosity curve of the debris at a par- 
ticular time step by computing the time at which the light 
was emitted from the debris and by finding the ionizing 
flux that was illuminating that location at the time the 
emission occurred, according to the light-travel time from 
the black hole to that particular region of the debris. 

Figure 4 shows three different Ha light curves from 
the debris confined to a plane (assuming = 500, 
^out — 40, 000) during the 60-day accretion phase of the 
5k simulation. Figure 4a shows the accretion luminosity 
on a logarithmic scale (solid curve), calculated from the 
accretion rate of the debris in the SPH simulation. The 
UV/X-ray luminosity curve is arbitrarily scaled and over- 
plotted on the top of the Ha curve for comparison. It is 
noticeable that the Ha light curve departs from the ac- 
cretion light curve at late times, though the departure ap- 
pears small in the logarithmic plot (used here due to the 
large dynamic range of the light curves). The same effect 
is more noticeable in the Figure 4b, where the accretion 
luminosity is proportional to t~^/^ and the Ha light curve 
is plotted on a linear scale. The Ha light curve roughly 
follows the shape of the incident UV/X-ray light curve at 
early times but decays faster at late times. The faster 
decay in the Ha light curve reflects the debris evolution 
in time: as the tail becomes more elongated, the incident 
photons travel a longer way to illuminate the debris. Con- 
sequently, the intensity of the illuminating light gets lower 
in the later stages of the tidal disruption event. The rel- 
ative decay rate of the Ha light curve with respect to the 
UV/X-ray light curve diminishes about 80 days after the 
accretion started. We find that this late relative rise in the 
Ha luminosity is due to the increase of number of parti- 
cles in the disk component. As particles diffuse from the 



high density tail to the lower density disk, in later stages 
of the simulation, their emission efficiency increases and 
they contribute a significant amount of Ha light to the 
light curve. To isolate the effect of the debris evolution in 
time from the evolution of the illuminating light curve, we 
calculate the Ha light curve in the case of constant illu- 
mination (Fig. 4c). Here, the relative departure of the Ha 
light curve from the UV/X-ray light curve can be inter- 
preted as a consequence of the expansion and redistribu- 
tion of the debris. The Ha hmiinosity appears to level off 
at late times because the debris disk begins to settle into 
a quasi-steady configuration. 

In summary, the observed Ha flux depends sensitively 
on the UV/X-ray light curve, on the distribution of mat- 
ter that makes up the inner portion of the debris, and 
on how quickly particles redistribute themselves in phase 
space. The main features of the Ha light curve are: an 
initial rise followed by a decline, with superposed fluctua- 
tions. The initial rise is a consequence of the propagation 
of the initial illumination front through the debris at the 
speed of light. The fluctuations are a result of the fluctu- 
ations in the accretion rate, which are caused, in turn, by 
the finite number of particles employed in the simulation. 
The decay rate of the Ha light curve is determined by the 
decay rate of the UV/X-ray light curve, debris expansion 
and redistribution rate. 

4.2. Line Profiles From the Debris and Their Variability 

We have computed sample line profiles emerging from 
the debris for the following choices of model parameters: 

1. Inclination angle of the plane of the debris, i. - We 
assumed that the observer is located on the positive 
z-axis, at a distance d ^ 00 at i = 30° to the z'- 
axis. Changing the inclination, changes the values 
of the projected velocities (i.e., the overall width 
of the line profile) but has very little effect on its 
shape otherwise. 

2. Inner and outer radius, and ^out- -The adopted 
inner radius of the debris is the inner boundary 
of the region from which Ha emission is expected 
to emerge. The choice of the inner radius de- 
pends on physical conditions in the debris, as we 
explain in §5.1. Here, we explore several cases 
with between 200 and 10,000. The outer radius 
^out = 40, 000 is naturally set by the dimensions of 
the system. 

3. Particle emissivity power-law index, (3. - As noted 
in §4.1, above, we find that = 2.4 in the debris 
tail and (3'^'^^^ = 0.06 in the disk. These indices de- 
scribe the emissivity per particle as a power-law and 
correspond to indices g*"*' and (f^^^ in the descrip- 
tion of the surface emissivity. The adopted values 
for emissivity indices significantly influence the pro- 
flle shapes. Higher values of (3 weigh the emissivity 
towards smaller radii where the projected velocity 
is higher. 

4. Velocity dispersion, a , and central source radius ^o- 
- The adopted value of the velocity dispersion for 
the profiles presented was a = 100 km s~^. A lower 
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limit on the velocity dispersion is set by the veloc- 
ity difference measured for pairs of adjacent parti- 
cles. This value represents the dispersion due to 
the finite number of resolving elements in the sim- 
ulation and equals 20 km for 90% of particles 
close to the debris plane. The constraint on the up- 
per limit of the velocity dispersion comes from the 
velocity dispersion due to small scale turbulence, 
c^turb < (l/2)(L/r)VKepier « 800 km s~\ where 
VKepier if' a local circular velocity, r is a distance 
from the massive black hole and L is the small- 
est dimension of the fluid (i.e. width or thickness 
of the tail). Since the real turbulence could be 
substantially smaller than the upper limit due to 
dissipation by internal shocks, we adopt a value of 
100 km s~^. Larger values of a produce wider pro- 
files and smooth-out sharp features. The central 
source radius is arbitrarily chosen to bo = 200. 
It implies a central source of finite dimensions such 
as a corona of ionized plasma, or a vertically ex- 
tended accretion flow in the innermost parts of an 
accretion disk. Its effect on the line profiles is rather 
small. 

In Figures 5 9 wo show sample lino profiles to illustrate 
how they evolve in time and how they are affected by the 
choice of model parameters and by the orientation of the 
observer. Figure 5 is a "trailed spectrogram" summariz- 
ing the temporal evolution of the line profiles from the two 
different SPH runs; it is a 2-dimensional map of the Ha 
emission as a function of projected velocity and time. Fig- 
ure 6 shows a different representation of the evolution of 
the line profile with time, which effectively comprises se- 
lected time slices from the trailed spectrogram. Figures 7 
and 8 show how the inner radius of the line-emitting re- 
gion and the azimuthal orientation of the observer affect 
the observed line profiles. Figure 9 shows the effect of 
the different values of velocity dispersion on the shape of 
the line profiles. The main properties and features of our 
results are as follows: 

Profile variability with time. - A property that is imme- 
diately obvious in the line sequence is the change 
of the profile shape with time (Figures 5 and 6). It 
is noticeable that the adopted low value of veloc- 
ity dispersion allows us to resolve individual par- 
ticles in the trailed spectrograms, orbiting around 
the black hole. The evolution of the line intensities 
in time roughly follows the behavior of the UV/X- 
ray luminosity but decays somewhat faster in time 
(see §4.1). The multi-peaked line profile is a conse- 
quence of the velocity field of the inner debris, which 
consists of the inner portion of the tidal tail that 
is falling towards the black hole (towards the ob- 
server) and debris that is rotating around the black 
hole after being scattered. The line profiles and 
their variability could be observationally important 
features of the debris just formed from tidal disrup- 
tion. The variable line profiles might be observed 
and recognized on the relatively short time scale of 
hours to days. 

Effect of the inner radius, . - The profiles become 
broader as the inner radius of the line-emitting re- 



gions decreases since higher-velocity gas resides at 
smaller radii (see Figure 7). The approximate full 
width at zero intensity of the profiles ranges from 
4, 500 km s"^ for = 10, 000 to 18, 000 km s"^ for 
Cin = 200. We find that line profiles change from 
the profiles dominated by the emission red-ward 
from the rest wavelength for < 1000 to nar- 
rower profiles dominated by the blue- ward emission 
from the tail for > 1000, since for large values 
of ^i„ , the high- velocity rotating gas in the vicin- 
ity of the black hole is excluded and the dominant 
contributions to the line profile come from the tidal 
tail. The intensity of the line also decreases with 
increasing inner radius, making the outer regions 
of the debris harder to observe. 

Effect of observer orientation, • - Because of the non- 
axisymmetric geometry and velocity field, the line 
profiles emitted by the debris, depend on the orien- 
tation of the tidal tail relative to the observer. In 
Figure 8 we show the effect of azimuthal orientation 
^0 of the debris, with respect to the observer. The 
values of (/-o arc 45°, 90°, 120°, 180°, 220° and 270°, 
as measured in a counterclockwise direction from 
positive x'-axis to the observer's line of sight. These 
can be compared with the profile corresponding to 
the same time in Figure 6 for (j)o = 0° . The position 
of the peaks in Figure 8 varies relative to the rest 
wavelength, since the relative direction of bulk mo- 
tion of the material depends on the observer's ori- 
entation. For example, it is possible to distinguish 
the emission from the tail for the range of azimuthal 
orientations 90° — 220°. The tail emission in these 
profiles appears as the most blueshifted peak, since 
these are the orientations for which different por- 
tions of the tail flow towards the observer. 

Effect of velocity dispersion, a . In Figure 9 we show the 
effect of four different values of the velocity disper- 
sion in calculation of the emission line profiles. A 
value of 20 km is the lower limit of velocity 
dispersion set by the discrete nature of the SPH 
simulation, while the upper limit of 800 km s~^ is 
determined by the small scale turbulence in the de- 
bris. The velocity dispersion of particles in the halo 
(about 1500 particles in the 20k run) is significantly 
higher and reaches 6000 km s~ . As the value of the 
velocity dispersion increases the profile features get 
smoothed out, until only a smooth, double-peaked 
profile is observed. 

Effect of illuminating light curve. - We have computed 
model profiles for several different X-ray illumina- 
tion light curves keeping all the other parameters 
fixed. We used (a) the light curve obtained from 
the accretion rate in the 5k simulation, (b) the light 
curve from the accretion rate as predicted by the- 
ory, i.e. oc t-^/^ (Rees 1988; Phinney 1989), and 
(c) a light curve that is constant in time (Figure 4) . 
We find that the line profile shapes do not depend 
sensitively on the shape of the light curve. This is 
a consequence of the centrally "weighted" emissiv- 
ity profile of the debris which causes the innermost 
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emission region to be the dominant contributor of 
the Ha hght. In the innermost region of the debris 
the dynamic range in hght-travel times is not large; 
therefore the illumination of the innermost emitting 
region is almost instantaneous. Over the very short 
light-crossing time of the central emitting region, 
the gradient in the UV/X-ray light curve is small 
and the illumination is nearly constant over this re- 
gion. The fast fluctuations in the illuminating light 
curve on the other hand arc smoothed out during 
reprocessing in the debris, and cannot be identified 
in the Ha light curve. 

In summary, we find that profiles are not significantly 
influenced by the shape of the illuminating light curve. 
The profile shapes are affected, however, by the inner ra- 
dius of the line-emitting region and redistribution of the 
debris. The inner radius can change with the advance 
or recession of the ionization front into the debris, which 
is controlled primarily by the density of the debris. The 
most notable effect of the inner radius is on the width of 
the profile. Since the physical conditions can change very 
rapidly during flares, this mechanism causes the line pro- 
flies to change on the light crossing time scale (minutes to 
hours) and evolve from wide multi-peaked to narrow and 
vice versa (the recombination time is negligible in com- 
parison to the light-crossing time of the debris, therefore 
particles in the debris respond effectively instantaneously 
to changes in the incident flux). The redistribution of 
the debris in phase space, on the other hand takes more 
time: ^24 hours for particles in the innermost part of the 
emitting region. This redistribution of the debris also may 
cause a transition from narrower to wider profiles, however 
it takes place gradually, on the time scale of days. 

5. DISCUSSION 

5.1. Physical Conditions in the Debris and Radiative 
Processes 

We have calculated the physical conditions and radia- 
tive processes in the debris using the photoionization code 
CLOUDY, version 94 (Ferland 1996). It is necessary to 
know the physical conditions in order to assess the validity 
of the assumptions made in our line profile calculations. 

Since the physical conditions in the tidal tail, disk and 
halo differ noticeably, these three regions naturally emerge 
as separate components of the tidal debris. The spheroidal 
halo is an oblate structure of particles scattered out of the 
plane of the debris. The disk is produced by the fiow 
of particles from the tail which turn around the black 
hole and form a higher density circular component con- 
centrated close to the plane of the debris. The disk shows 
a smooth transition to the halo in terms of density and 
physical parameters (see Table 2). The physical differences 
among the components arise as a consequence of the num- 
ber density, which is several orders of magnitude higher in 
the tail. The particle density of the halo is uncertain, due 
to the numerical scatter, and consequently its luminosity 
contribution is also subject to uncertainty. We neverthe- 
less have calculated and presented physical properties of 
the spheroidal component, as given by the SPH simula- 
tion, and we discuss the implications of its presence for 
the line profiles and the observability in the next section. 



The temperature of the debris tail reaches 3 x 10^ K in 
the hottest parts of the tail (i.e. the illuminated face of 
the tail) and equals 5000 K on average in the partially ion- 
ized and neutral parts of the tail. The temperature in the 
disk ranges from 8000 K in the inner region to 1 x 10^ K 
at the outer rim of the disk, with the average ~ 10^ K. 
The mean temperature of the halo is 1 x 10* K and ranges 
from 3 X 10^ K in the central parts of halo, to 1 x lO^'K 
in the outer halo. The ionization parameter is calculated 
for all components as U = ^n/nn c, where '^h is the flux 
of ionizing photons, hh is the total hydrogen density, and 
c is the speed of light. The ionization parameter in the 
tail ranges from 0.1 in the parts of the tail closest to the 
source of ionization, to 10^'"' on the far side of the tail. 
The ionization parameter in the disk is almost constant 
throughout the disk with a value of about 20. The ioniza- 
tion parameter in the halo is higher with an average value 
of about 27. As a consequence, the disk and the halo arc 
in a much higher state of ionization relative to the tail. 
The disk shows a wide range of hydrogen ionization frac- 
tions over its radius, with the strongest ionization at the 
outer rim of the disk where density is lowest. In the halo 
the fraction of neutral hydrogen atoms ranges from 10""^" 
to 10~^, with an average value close to 10"''. In contrast, 
the tail is mostly neutral, with only one hydrogen ion per 
10^ hydrogen nuclei. The properties of the tidal disruption 
regions are summarized in Table 2. 

The ionization state of the debris as well as its abun- 
dance and density determine the radiative processes dom- 
inant in the debris. The implication of the assumed solar 
metalicity is that emission lines from metals play an im- 
portant role in the cooling of the debris. We focus our 
attention on processes relevant to the final Ha luminos- 
ity. Two processes contribute to the Ha luminosity of 
the debris: recombination and collisional excitation. The 
dominant recombination channel is the recombination of 
a photoclectron to the n = 3 level followed by a decay to 
the n = 2 level via emission of an Ha photon. Emission by 
collisional excitation occurs when hydrogen atoms in the 
n = 1 and n = 2 levels are promoted to the n = 3 level 
by collisions with energetic photoclcctrons and then de- 
excite radiatively. The relative contribution of Ha emis- 
sion through recombination relative to collisions is 

-^Ha,rec ^rec I^HII ^rec 

-^Ha,coll acoll nHI Vcoll 

where a^ec and acoii are the respective coefficients for re- 
combination and collision processes which lead to emission 
of an Ha photon, nnii and nni are the number density of 
hydrogen ions and atoms, and Vrec and Vcoii arc the parts 
of the debris volume which gives rise to emission through 
recombination and collisional de-excitation, respectively. 

The recombination coefficient for an electron to re- 
combine from the continuum to the n=3 level can be 
writcn as ac = X]^3 '^n fo'' a given temperature. The 
adopted recombination coefficients for the tail and halo are 
a*°j? = 3.4 X 10-13 cm3 s'^, a^f ~ 1.8 x lO'^* cm^ s'S 
and a^f^° = 1.8 x IQ-^^ g-i (Qsterbrock 1989). Sim- 
ilarly, it is possible to estimate the rate of collisions which 
will promote an electron to n=3 state. The high opti- 
cal depths in the Lyman series cause a significant elec- 
tron population in the n = 2 level. The photoelectron 



Tidal Disruption Signature 



9 



energy distribution for low-energy levels of hydrogen dif- 
fers from the Boltzmann distribution and the number den- 
sity is not sufficiently high for collisions to lead to ther- 
mal equilibrium (Eq. 10.57 of Krolik 1999). In partic- 
ular, thcrmalizcd energy levels will be n>l for the tail, 
n>2 for the halo, and n>3 for the disk, and the rela- 
tive populations of these levels will be in agreement with 
their statistical weights. It is possible then to estimate 
the coUisional excitation coefficients for transitions from 
the n==l and n=2 levels to the thcrmalizcd n=3 level 
(Eq. 3.21 of Osterbrock 1989). This approximately yields 
coUisional excitation rates for the tail, disk and halo of 
a*«i| = 1.2 X 10-" cm'^ s^^, af^ff = 1.1 x IQ-^ cm^ s'^, 
and a^^/;° = 9.0 x 10"" cm-"* s-\ 

With the values of the rate coefficients, the ionization 
state of the tail, disk and halo {riHii/nH ~ 10"'^ in the 
tail and lO'' in the halo and disk on average), and taking 
into account the ratio of emitting volumes for two pro- 
cesses, it is possible to calculate the ratio of respective rel- 
ative contributions from recombination and collisions from 
equation (9) to be about 1 x 10"^, 6 x 10^, and 2 x lO'^. 
These vahies arc not unexpected; coUisional excitation is 
the dominant mechanism for production of Ha photons in 
the tail, where the density is highest, while recombination 
is dominant in the disk and halo. 

Some of the Ha photons created by the two processes are 
destroyed on the way out of the debris. The mechanisms 
that impact the total output in Ha luminosity are the ab- 
sorption in the Ha transition, and electron scattering, with 
respective optical depths THa and Tg. The optical depth 
"^tia very high for the Ha photons traveling across the 
tail, along the lines of incidence of illuminating photons 
and typically has a value of order 10^. Electron scatter- 
ing in the tail has an average optical depth of r*'**' « 10. 
Consequently, the majority of Ha photons are created in 
the thin ionized and partially ionized layers of the tail (we 
refer to them as the tail "skin"). From here a certain 
fraction of Ha photons escapes the debris and reaches the 
observer. In the disk rf}^ w 4 x 10® and rf^"^ w 8. In 
the early stages of accretion, while the luminosity is still 
super-Eddington. the halo is optically thin to Ha pho- 
tons (t]^'Jj'° ft: 2 X 10~^) and the electron scattering optical 
depth is about r^"'° » 5. In the late stages of the debris 
evolution r'jf^" » 1 x 10** and Tg stays approximately the 
same in the halo. The high values of r of the halo mate- 
rial destroy the fraction of the Ha photons emitted by the 
part of the debris embedded in the halo. The evolution of 
the spherical halo would further help the process of fad- 
ing of Ha emission with time along with the decreasing 
accretion luminosity. For the purposes of radiative trans- 
fer calculations we have modeled the tail, the disk and the 
halo as three separate components and therefore we do not 
account for the secondary absorption and electron scatter- 
ing by the halo, of Ha photons created in the tail. Both 
opacity mechanisms cause the destruction of Ha photons, 
if the spherical halo is fully evolved, and possibly wipe out 
the Ha emission line. If the same process operates in the 
centers of galaxies which are tidal disruption candidates, 
it should be possible to observe the disappearance of the 
broad Ha emission line on the scale of months. 



5.2. Spheroidal Halo: Implications for Line Profiles and 

Luminosity 

In our predictions of the observational signature of the 
post-disruption debris we have not included any contribu- 
tions from the spheroidal halo because we doubt that this 
structure, as predicted by our SPH simulation is real. In 
this section we discuss this issue further and justify our 
approach. We also describe qualitatively the effect of such 
a halo on the observational appearance of the debris. 

The halo is very likely produced artificially because a 
number of particles is scattered from the disk every time 
the Uow of particles turns around the black hole and inter- 
sects with the inflowing stream of particles. The munber 
of particles contributed to the halo by numerical scatter 
is proportional to ViV, where N the total munber of par- 
ticles in the run, and multiplied by a factor Trun/Tdyn, as 
described in §2. The dynamical time scale of the disk is 
about 7 days for the outermost particles, which implies 
that the number of scattered particles is of the order of 
1000. This is an upper limit of the uncertainty, where to- 
tal number of particles in the halo, in the 20k run, reaches 
number of 1500 particles in the later stages of the debris 
evolution. 

It has been pointed out by several authors that a 
spheroidal halo (Ulmer, Paczynski & Goodman 1998; Locb 
& Ulmer 1997; Cannizzo et al. 1990; Alexander & Livio 
2001) with a central radiatively supported torus (Rees 
1988; Loeb & Ulmer 1997; Cannizzo et al. 1990; Evans & 
Kochanek 1989) may form around the black hole as a con- 
sequence of a self-compression of the flow, which leads to 
the ejection of the debris with high velocities, perpendic- 
ular to the plane. Howcn'cr, none of the above mentioned 
structures were resolvcid beyond doubt in the simulations. 
A halo-like structure is expected to be in hydrostatic equi- 
librium, and its size should be determined by the radia- 
tion pressure from the central source. Ulmer, Paczynski 
& Goodman (1998) and Locb & Ulmer (1997) argue that 
the luminosity incident upon the halo can rarely be exactly 
tuned for the structure to rest in equilibrium. They further 
calculate that for a super-Eddington luminosity the halo 
will expand and cool down on the time scale of months 
to years, until it becomes gravitationally unbound and is 
blown away. For a sub-Eddington luminosity, on the other 
hand, the halo is expected to collapse due to insufficient 
support from radiation pressure. 

Based on its physical conditions it is obvious that, if 
real, the halo observed in our simulations would make the 
dominant contribution to the line emission from the debris. 
In the late stages of the debris evolution, our simulations 
show the halo becoming optically thick to Ha photons. 
Moreover, the long diffusion time of photons in the halo 
{tdifj ^ Te ^/c ~ 3.5 hours, where r the is inner radius of 
the halo) may smear out the time variability of the emis- 
sion line profiles, which is one of the main signatures of 
the disruption event. The high value of the velocity dis- 
persion of the particles in the halo would smear out the 
line profiles and very likc;ly make them unobservable. It is 
therefore necessary to address in future studies what frac- 
tion of the halo (if any) is really present and what fraction 
is contributed by the numerical scatter. 

If a halo manages to form, survives the super-Eddington 
phase, and achieves hydrostatic equilibrium it will be 
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transparent to Ha photons only during the super- 
Eddington phase. At later times, we expect that it will 
become optically thick and will be the primary source of 
illumination of the outer debris. Moreover, it will fade on 
a time scale of months to years. In this picture, which is 
consistent with the predictions of Loeb & Ulmer (1997), 
the shape of the observed line profiles should be similar 
to the profiles computed here for large values of the in- 
ner radius (see Figure 7). Whether the profile (noticeable 
here for ^ = 10, 000 at the wavelength A = 6500l) will 
appear as blueshifted or redshifted would depend on the 
orientation of the debris tail with respect to the observer. 

5.3. Observability of Emission Lines and Their 
Uniqueness as a Tidal Disruption Signature 

The CLOUDY calculation predicts a time-average Ha 
luminosity from the tidal debris of about 10^^, 10^^, and 
6 X 10^^ erg for the tail, disk and halo, respectively. In 
the earlier stages of the disruption event when the UV/X- 
ray luminosity is super-Eddington, the Ha luminosity is 
expected to be up to 80 times higher than its average value 
and comparable to that of tidal disruption candidates ob- 
served in the local universe (see Figure 4). The examples 
are NGC 4450 (at 16.8 Mpc) with an Ha luminosity of 
LHa=1.8xl039 erg s"^ (Ho et al. 2000) and NGC 1097 (at 
22 Mpc) with Ljy„=7.7xl0^^ erg s"^ (Storchi-Bergmann 
et al. 1995). Thus the emission-line signature of a tidal 
disruption event should be detectable at least out to the 
distance of the Virgo cluster. In practice, however, the de- 
tection of such emission lines from low luminosity sources 
may be complicated by the weak contrast relative to the 
underlying stellar continuum. 

Based on observational constraints from known tidal dis- 
ruption candidates, it should be possible to detect some 
variable properties of the Une profiles and light curves pre- 
dicted here. One of the first observable effects of tidal 
disruption should be UV/X-ray flash accompanied by a 
decaying light curve, mirrored in the delayed response of 
the Ha light curve of the debris, with some scatter. The 
line profile intensities are expected to decay accordingly 
in time. Other effects to look for are the change in the 
number of peaks in the line profile, relative fluctuations in 
intensity of the peaks as well as their shift in wavelength. 

The temporal variability of the Ha emission line profiles 
from the post-disruption debris is one of the important in- 
dicators of a tidal disruption event. In order to capture 
the rapid profile variability, due to the variable illumina- 
tion, the exposure time should be comparable to the light 
crossing time of the innermost regions of the line-emitting 
debris, which has the fastest and strongest response to the 
ionizing radiation. Longer exposures are expected to cap- 
ture the average shape of the rapidly varying line profiles. 
The light- crossing time of the innermost regions of the de- 
bris is about 8 ^2 Mq minutes (where ^2 = ^/lOO), while 
the exposure times are typically about 30-60 minutes (for 
galaxies at the distance of the Virgo cluster, for example). 
Thus, if an event is caught early in its evolution and the 
light-crossing time is relatively long (i.e., M > 10^ M0), 
there is a chance of detecting variability caused by chang- 
ing illumination over the course of one to a few nights. On 
longer time scales, variability is caused by changes in the 
structure of the debris. In the presence of the spheroidal 



halo, the variability of the lines may be modified by the 
long diffusion time scale of photons through the halo. The 
component of the tidal tail outside the halo will then still 
respond to the variability but on the time scale set by the 
light reprocessed by the halo. 

In view of the predictions from our profile calculation, 
an important question is whether variable multi-peaked 
line profiles can originate in some other physical scenario 
or can be regarded as the unique signature of tidal dis- 
ruption. Multi-peaked emission lines are likely to be the 
signature of an inhomogeneity in the phase-space distri- 
bution of the emitting material. Due to the asymmetry 
of the emitting region, the direction of the observer has 
a large infiuence on the observed line profile, making it 
difficult, if not impossible, to infer the exact emission ge- 
ometry from a particular multi-peaked profile. Neverthe- 
less the variability pattern of the line profiles can serve 
as a general indicator of a tidal disruption event. Addi- 
tional observational indicators can be used in conjunction 
to diagnose a tidal disruption event, for example, a sharp 
X-ray/UV flash preceding the appearance of the emission 
lines, an emission line spectrum indicative of a hard ioniz- 
ing contimmm (i.e., the presence of ionic species with high 
ionization potentials), and the characteristic decay of the 
emission-line flux on time scales of weeks to months after 
the event (the flux of different emission lines is expected to 
decay at different rates; see, for example, Eracleous, Livio, 
& Binette 1995). 

5.4. Approximations in the Calculation 

The finite resolution of the simulation could introduce 

uncertainties in the accretion light curve caused by the 
discretized accretion of the stellar material onto the black 
hole. Another possible effect is that the finite number of 
particles does not completely reveal the morphology of the 
debris and some of its components may stay hidden. For 
example, when the majority of the particles are confined 
to the tail it is hard to say if there is an accretion disk 
forming around the black hole out of a small number of 
particles. Consequentially, there is a concern that profiles 
calculated for low values of the inner radius, i.e. ^i„=200, 
may represent the contribution from a small number of 
particles in the central region and therefore dominated by 
small-number noise. 

A possible source of error is the assumption of a thin, 
flat structure (i.e., confined to a plane). In the early stages 
of the evolution of the debris, the majority of the particles 
are still in the tidal tail, located in a plane, which makes 
the assumption valid. During the evolution of the debris 
the number of particles that orbit around the black hole 
in an almost spherical distribution increases (Cannizzo et 
al. 1990; Loeb & Ulmer 1997; Ulmer, Paczynski & Good- 
man 1998; Ulmer 1999; Menou & Quataert 2001). At that 
point it is possible to distinguish three structural compo- 
nents of the debris: the relatively planar tidal tail and 
the disk and the spheroidal halo. The halo is made up of 
particles scattered from the tail by shocks during the peri- 
centric approach of the debris or during the intersection 
of the tail with itself (Kochanek 1994; Lee, Kang, & Ryu 
1996; Kim et al. 1999; Ayal et al. 2000) and some fraction 
of particles contributed by the numerical noise. Since the 
tidal tail includes the majority of the particles, and most 
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of the halo mass is concentrated close to the equatorial 
plane, even in our last frame, the assumption of a thin 
disk is still reasonable. If, however, the mass in the spher- 
ical halo increases at very late times, the assumption of a 
planar geometry needs to be reconsidered. 

We also adopt the weak field approximation in our cal- 
culations, which is a fairly small source of error (of order 
1% and less). This is valid since we adopt = 100 as 
the innermost radius of the line-emitting debris. The few 
particles interior to this radius would not contribute to the 
Ha emission because their close proximity to the source of 
ionizing radiation would make them fully ionized. In the 
case of a physical scenario where emission of Ha is not 
possible because of a highly ionized debris, the same pro- 
file formalism can bo used to calculate the emission from 
other lines emitted under these conditions. 

6. CONCLUSIONS 

We modeled the emission-line luminosity and profile 
from the debris released by the tidal disruption of a star 
by a black hole in the early phase of evolution. Our model 
predicts prompt optical evolution of post-disruption de- 
bris and profile shapes different from circular and elliptical 
disk model profiles. Since line profiles observed so far in 
LINERS look more disk-like and evolve slowly, the obser- 
vations are likely to have caught the event at late times 
(> 6 months after the initial disruption), after the debris 
has settled into a quasi-stable configuration. 

The line profiles can take a variety of shapes for dif- 
ferent orientations of the debris tail relative to the ob- 
server. Due to the very diverse morphology of the debris, 
it is almost impossible to uniquely match the multi-peaked 
profile with the exact emission geometry. Nevertheless, 
the profile widths and shifts are strongly indicative of the 
velocity distribution and the location of matter emitting 
the bulk of the Ha light. Profile shapes do not depend 
sensitively on the shape of the light curve of the X-rays 
illuminating the debris. They strongly depend on the dis- 
tance of the emitting material from the central ionizing 
source, which is a consequence of the finite propagation 
time of the ionization front and the redistribution of the 
debris in phase space. It may be possible to distinguish 



between the two effects observationally, based on their dif- 
ferent characteristic time scales. The onset of the optically 
thick spheroidal halo should cause the disappearance of the 
broad Ha emission line on the time scale of months, and 
give rise to the emission of narrower, strong, blucshiftcd 
or redshifted emission line, arising from the portion of the 
tidal tail unobscured by the halo. 

If X-ray flares and the predicted variable profiles could 
be observed from the same object they could be used to 
identify the tidal disruption event in its early phase. The 
X-ray flares can be promptly detected by all-sky synoptic 
X-ray surveys and high energy burst alert missions such 
as Swift. The evolution of the tidal event may then be 
followed with optical telescopes from the ground on longer 
time scales and give an insight in the next stage of de- 
velopment of the debris. Thus, simulations of the tidal 
disruption process on longer time scales (of order several 
months to a few years) are sorely needed. Calculations 
of the long-term evolution of a tidal disruption event can 
predict the type of structure that the debris finally settles 
into and whether its emission-line signature resembles the 
transient double-peaked lines observed in LINERs. This 
study would provide an important insight into the evolu- 
tion of LINERs. 

Finally, the observed rate of tidally disrupted solar type 
stars can constrain the rate of captured compact objects 
(which are important gravitational wave sources), and the 
capture rate of main sequence stars in our Galaxy, which 
are expected to emit the peak of the gravitational radia- 
tion in the LISA frequency band and can be detected in 
the local universe. 
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Table 1 

Paraaieters of SPIl Rrxs 







Total 


Start of 


Emission-Line 


Time 


Particle 


SPH 


Total 


Duration 


Accretion 


Evolution 


Step 


Mass 


Run 


Particles 


(days) 


(days) 


(days) 


(hours) 


(g) 


5k 


5000 


94 


34 


60 


4.11 


3.978x10^9 


20k 


20000 


53 


47 


6 


0.746 


9.945x10^8 



Table 2 

Physical Properties of the Debris 





Nh 




riHi/nHtot 


nHii/riHtot 


2^ b 


U ^ 


r d 


Debris Region 


(cm^2) 


(cm-3) 






(K) 




(erg s-i) 


Tail 


1027-1031 


10l4_10l5 


~ 1 


- 10-3 


5 X 103 


10-5-0.1 


1 X 1036 


Disk 


1021-1025 


1011-1012 


10-8-0.3 


0.7-1 


~ 1 X 10^ 


20 


1 X 1037 


Halo 


1025 


10^-1012 


-10-7 


~ 1 


1 X 10^ 


27 


6 X 1038 



'^The column density: low and high values correspond to directions orthogonally and radially through 
the debris component. 

''The average value of temperature over radius. 

'^Ionization parameter: range of values in the tail and average in the disk and halo. 

'^Tlie Ho luminosities from the debris components as calculated for the time-average value of incident 
luminosity Lace = 1-5 x lO'^^eig s~^. 
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Fig. 1. — Maps showing the positions of the SPH particles from the 5k simulation after the second pericentric passage, projected in the 
x'y'-plane at two different times. Particles gravitationally bound to the black hole arc colored black, while the unbound particles arc colored 
grey. The dashed line represents the initial trajectory of the star before disruption and the trajectory of the center of mass of the debris 
after disruption. The maximum particle velocities are of order 10^'^ c. Upper Panel: Particle map at the start of the accretion phase, 34 days 
after the disruption occurred. Inset: Particles in the inner region of the debris, orbiting close to the black hole. Lower Panel: Particle map at 
the end of the simulation, 94 days after the disruption. The tidal tail can be clearly separated into particles that are unbound and about to 
escape the black hole and particles that arc returning towards the black hole. The inner region of the debris consists of returning particles 
from the inner tail that have been scattered and form a rotating structure around the black hole. 
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Fig. 3. — The quantum efficiency for partictes in the dislc (upper) and the taii (tower) as a function of distance from the black hole. The 
data points are the results of CLOUDY photoionizations calculations performed for different regions in the tail and the disk. The solid line 
is the best fit power law, with power law index, q*''''=2.41 and g''*^'^=0.06, respectively. 
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Fig. 4. — The Ha light curves (data points) resulting from the re-processing of three different illuiniiiatioii light curves. The solid line 
shows a scaled version of the UV/X-ray continuum light curve that illuminates the debris, which in a) follows the SPH accretion rate; b) 
decays with time as t~^/^; c) is constant. Note that the vertical axis is calibrated logarithmically in (a), while in (b) and (c) they have the 
same linear scale. The Ha light curves rise initially as the illumination front propagates through the debris and then decay faster than the 
UV/X-ray light curves. The details of the calculation are described in §4.1. 
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Velocity (km s ') Velocity {km s ^) 

-10^ 10"* 2x10* -10* 10* 2x10* 




Wovelength (A) Wovelengtli (A) 

Fig. 5. — Trailed spectrogram of the simulated Ha emission-line profiles from the 5k simulation spanning 60 days (left) and from the 20k 
simulation spanning 6 days (right). This is effectively a 2-dimensional intensity map versus projected velocity of the emitting material and 
time. Darker shades correspond to higher intensity. The scale on the right represents time since the tidal disruption event. 
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Fig. 6. — Sequence of Ha profiles emitted from the region ^ G (500,40 000) over period of 6 days (20k rim). The relative time from 
the beginning of the accretion phase onto the black hole is marked next to each profile. The accretion phase begins 47 days after the tidal 
disruption. The inclination of the debris plane and the velocity shear are as marked on the figure. 
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Fig. 7. — Ha emission line profiles simulated for seven different values of inner radius (§in)i as determined by the propagation of the 
ionization front through the debris. The relative time for profile frames is S'* 6^ 0™. The intensity of the profile calculated for = 10, 000 
is multiplied by the factor of 100. The inclination and velocity shear are as marked on the top of the figure. 



Tidal Disruption Signature 



19 




Fig. 8. — Ha emission line profiles simulated for six different azimuthal orientations of the debris with respect to the observer, as represented 
by 00. See profile in Figure 6 with the time label S'' 22*^ 30"" for orientation iPq = 0° . The size of the emitting region, inclination and velocity 
shear are as marked on the top of the figure. 
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Fig. 9. Hq emission line profiles simulated for four different values of velocity dispersion. A velocity dispersion of 100 km s^^ is the 
equivalent of an instrumental resolution of 2 A. The relative time for the profile frames is 6'* 6^ 0™. Size of the emitting region and inclination 
are as marked on the top of the figure. 



